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Abstract 

Using the first-order approximating solutions to the Einstein-Maxwell- 
Klein-Gordon system of equations for a complex scalar field minimally 
coupled to a spherically symmetric spacetime, we study the feedback 
of gravity and electric field on the charged scalar source. Within a per- 
turbative approach, we compute, in the radiation zone, the transition 
amplitudes and the coherent source-field regeneration rate. 
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Boson stars, generally seen as gravitationally bound, both globally U(l) 
and spherically symmetric, compact equilibrium configurations of cold com- 
plex scalar fields, were discovered theoretically, more than 30 years ago, by 
Kaup [1] and Ruffini and Bonazzola [2]. Later, a major interest has been 
focused on macroscopic stable boson stars, once they arose as promising can- 
didates for non-baryonic dark matter in the universe [3]. Moreover, according 
to current observations, it has been suggested that a supermassive boson or 
soliton star could be at the center of our galaxy [4]. These assumptions 
raise the question whether such configurations, if they exist, are dynamically 
stable [5] and the settling down of a boson star to a stable state has been 
investigated via numerical calculations [6]. 

In our approach, we term by a boson star the charged scalar nebula which 
finds itself in one of the spherically symmetric positive-frequency modes of 
radial wave number k, with k 2 — > 0+, and pulsation Uk = [m^ + k 2 } 1 ^ 2 . Such 
a configuration is obviously unstable and the instabilities is expected to lead 
to the formation of boson star from a initially smooth state [7]. Not very 
much is specifically known in this direction [8], although the reversed stability 
to instability passage has been extensively investigated [5, 6, 9]. 

In what it concerns an analytical approach, exactly solvable models for 
boson stars with large selfinteraction [10] and for boson-fermion stars [11] 
have been worked out only in low dimensional gravity. In four dimensions, 
the bosonic or the mixed fermion-bosonic fields interacting via gravity have 
been investigated mainly by numerical calculations [1-6,8,9,12]. Therefore, 
the analytical solutions of the coupled field equations could be of interest 
for a better understanding of different stellar configurations as well as for a 
numerical-functional combined iterative treatment which describes the dy- 
namics of charged boson stars formation. 

In a previous letter [7], we have studied the 5*0(3, 1) x [/(l)-gauged mini- 
mally coupled charged spinless field to a spherically symmetric spacetime and 
analytically obtained the first-order approximating solutions to the system 
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of Klein-Gordon-Maxwell-Einstein equations. These allow us to go further 
and to study now the feedback of gravity and electric field on the charged 
scalar source, computing perturbatively, in the long range approximation, 
the correspondingly induced transition amplitudes. 

Let us consider the spherically symmetric configuration described by the 
metric 

ds 2 = e 2 f {drf + r 2 (d9 2 + sin 2 9 d<p 2 ) - e 2h (dt) 2 , (1) 

where / and h are functions of r and t, and define the pseudo-orthonormal 
tetradic frame {e a } a= jj with the corresponding dual orthonormal base 



UJ 



°J dr , 



rdd, u 3 = r sin 9 dip , uj a = e h dt 



(2) 



The 5*0(3, 1) x [/(l)-gauge invariant Lagrangean density for a charged 
boson of mass m , coupled to the electromagnetic field, 



C = V ab M;b + m 2 00 + -F ab F ab , 



(3) 



lead to the following system of Klein-Gordon-Maxwell-Einstein equations 

□0 - m\<\) = 2ieA c (p\ c + e 2 A c A c (p , and its h.c. , (4) 

(5) 
(6) 



F ab . h = 



te rj 



ah 



(0| 6 - ieA h (f)) - (cf)\ b + ieA h (f)j <p 



Gab = K (0;a0;6 + 4>;b4>;a + F ac F b C - r] ab C^j 

The components of the Einstein tensor, G a b, have the explicit form 



G22 



G33 = 



1 - e~ 2 f 



h\ll + (/^|l) 2 - /|44 + (f\if 



G44 

G14 



= -2 



-2/ 



+ 



+ fyi 

r 

0-2/ 



r 



-2/ 



+ 



2 — /|4 

r 



(7) 
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where ( • )k = e <j( ' )• Also, in (3), <f)- a = <p\ a — ieA a <f) and the Maxwell tensor 
F a b = Ab :a — A a: b is expressed in terms of the Levi-Civita covariant derivatives 
of the four-potential A a , i.e. A a:0 = A a \ — A c T c ab . 

In our previous letter [7], we have assumed that the charged scalar field 
is the main source of both the electromagnetic and gravitational fields and 
neglected (at large distances and within the framework of a first-order approx- 
imation) the feedbacks of gravity and electromagnetism. The corresponding 
equation of motion has simply become the one of an spherically symmetric 
state on a Minkowskian background, i.e. 

2 

0,rr H — <f>,r ~ 4>itt — ™>o = 0, and its h.c. , (8) 
r 

with the positive-frequency mode solutions 

r r 

where Uk = [k 2 + m 2 ] ] 1 ^ 2 . Moreover, by neglecting the gravity feedback on 
the Maxwell sector also, as it effectively involves second-order contributions 
of the charged scalar 0, the Lorentz condition and the Maxwell equations (5) 
are satisfied, in the minimally symmetric ansatz A\ = Ai(r,t), A4 = A±(r,t), 
4> = 0(r, t), by the solution(s) 



A 1 = ek\Aff 



r _ . \M\ 2 



A 4 = 2euj k \Af\ 2 log— + 2ek^^t , (10) 

r r 

which correspond to the electric field 

E = Fu ~ A^ r -A u = 2eu; k l —^-2ek l -^-t (11) 

In the particular case k — 0, where the star is just above the passage to 
the possible stable excited states, its mode pulsation uj = m being located 



4 



at the accumulation point of the eigenfrequencies of an excited boson star, 
we have found for the linearized Einstein field equations 



h,r 9 f — K 



9 2 2 M 

2e m. 



h )Tr -\- {h^r f ir ) 



\M\ 



A 2 2 IA/T , r 



|AT| 



log — 



2 2 



the solutions 



2m, 



IM 2 , |A/"| 



r 2 



+ 



+ 2e 2 mg 



+ 4e 2 ml 



log 



(12) 



f(r) = 



C x b r 
+ a log — 



Ci 
r 



where 



/i(r) = 

b = k \J\f\ 2 , a = m 2 Q b , and |7V| 2 



1 

2^2' 



(13) 



(14) 



that have led to analytical expressions for total charge, particle number, 
radius and mass of the analyzed configuration [7]. 

In the followings, we are going further and analyze the feedback of gravity 
and electric field, respectively expressed by the metric functions (13) and 
the four-potential (10), on the Klein-Gordon equations (4). Within a first- 
order perturbative approach, we write down the wave function describing the 
charged scalar field as 



$(r,e,tp,t) = </>(r,t)+x(r,e,(p,t), 



(15) 



where \x\ <C \<f)\ and is the Minkowskian background solution (9), with 
|jV| = 1/V2e 2 . Consequently, the Klein-Gordon equations (4), with k — 
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implying uo = m , turn into 



d 2 x | 2 dx 
dr 2 r dr 
2Ci j. H£i «9x 
<9r 



1 _ 2c x d 2 x 2 



l-log^) 



+ e 



— 2im§e r log 



•9% 
r J dt 



£1 \ 



er 



er 



1 + log 




log 



r 



(16) 



and its h.c, and we shall focus on the long range behaviour, r < r < 00, 
where the above equation gets the simpler form 



2M 

e r 



d 2 x 2 dx 
dr 2 r dr 



I ; 2M <9 2 x 



+ - A X - e 



dt 2 



m ( r \ dx o, 2 ( r \ 
+ 2 imo er logy - + m log (-j x 



V2 



^log 2 - e 



r 



-irngt 



(17) 



The constant Ci has been replaced by the total gravitational mass of the 
Bose star, given by the Tolman's relation [13] 



M 



J T 44 e f+h Anr 2 dr , 



(18) 



and explicitly reading [7] 



M 



+ 



2tt 
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1 — a 



3 + log - — PolyGamma ^0, — 
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(19) 



Since, in the radiation zone, one can consider 



,±- 



2M 
1± + 



2M 



n>2 
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the equation (17) does actually become 



□X = 



2M { I d 



r 2 dr 



— 2im log 



dr 



r \ dx 



+ 



dt 2 J 

J?L e -irn t log 2 ( ' 



(20) 



To/ dt ^/2er Vr , 

where □ is the usual d'Alembertian on R 4 , and therefore, once it achieved 
the standard form, 

D X = V X + J, (21) 
employed in Perturbation Theory, it points out the operators 



V(r,t) = 
J(r,t) = 



2M { 1 d 



,2 



dr 



d 2 ) n . , fr\d 

+ -a- ^oiogy - 



(22) 
(23) 



\/2er \r , 

describing the (perturbed) effective potential and current. The massless-like 
0-th order equation 

d 2 h 



1 d f 2 dh\ 1 ~ 



dr 



8t 2 







(24) 



provides the complete orthonormal set of positive-frequency modes (in terms 
of spherical Hankel functions) 



(25) 



which enables us to compute the first-order transition amplitudes between 
the initial and final states as: 



A"jy = J K, llm ,(x) (VK lm (x))r 2 drdndt, 



(26) 



where 



Vh ulm (x) = -2 



M 



+ 2m OJ log 



r 



Ki m {x) , (27) 



with (x) = (r,9,if,t). In (26), we identify the transitions of the charged 
boson in the presence of gravity 



•4sO(3,l) 



A 11 



= -Muj z 
2 7 



H\%{Jr)H)'Uujr)dr 



r->0 



1+ 



(1) 
l+l 



Mu , for I = 0, 1 

7T 

Mi(i + r 



— H l+ Aur)H l+k 



r^O r 



Muo 



7T 



(28) 



and of the electric field potential 



*4{/(i) 



-m uj I H^ l ( UJ ' r ) 



r log ( — 

Vr 



(lot) dr 



UJ 



(29) 



where 7 is the Euler's constant 7 ~ 0.577216 and the well-known 5-Dirac 
multiplier 2tt5(uj' — uo) has been factored out. 

The corresponding transition rate for each of the above amplitudes is 
respectively given by 



(30) 



Finally, let us turn to the current operator (23) in order to study the 
possibility of spontaneous creation of charged bosons in the presence of the 
electric field potential A4. These process is described by the transition am- 
plitude 

J — [ Vuh* Ljlm (x)J(r,t)r 2 drdttdt 

.2 



A 



= — 2tv5(lu — m ) \/lix<jj — - 
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m5 

e Jr^O 



VtH[%(u;r) log 2 



(31) 



Integrating by parts with a cut-off in z 
becomes 

Aj = 2-k5(uo 

where the integral 

^ dz 



r/ro, the previous expression 

(32) 



mo) —Z, 



euo 



s: 



-ibz 



Z 



log(z) 



ir 2 7 2 i 

24 + Y + 2 77r + 



3^4 



in log(fo) 

7 H 

2 2 

4 



log(6) 
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(33) 



with b = ujt , contains the generalized hypergeometric functions 
p F q [{a±, . . . , a p } , {b±, . . . , b q } , which are finite for all finite arguments and 
p < q [14]. In the above calculations we have assumed that uj = 2n/T is the 
natural dual of the "repeating time-interval" such that 

rT r dt du 



lim 

T^oo Jo 



R+ 27T 



= E 



states 



1, 



and then, 



jV+{m Q -J) = 2tc5(lu 



mo) 



'2m 2 Q 
euo 



IXI' 



(34) 



does concretely represent the coherent source-field regeneration rate. 
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